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Various scaling relations of the entanglement entropy are reviewed. Based on 
the scaling, I would like to point out similarity of mathematical formulation 
among recent topics in wide research area. In particular, the scaling plays 
crucial roles on identifying a quantum system with a physically different clas- 
sical system. Close connection between the scaling and hyperbolic geometry 
and contrast between bulk/edge correspondence and compactification for the 
identification are also addressed. 
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1. Introduction 

The entanglement entropy is one of the most fundamental concepts in quan- 
tum information. Recently, its efficiency as a tool to see through underlying 
physical principles in our targets is also recognized in statistical physics, 
condensed matter physics, and string theory. This wide applicability comes 
from universality of the entropy irrespective of their details. Roughly speak- 
ing, the entropy represents a logarithm of a correlation function. Thus, the 
entropy directly picks up critical exponents. In addition, the entropy can de- 
tect topological structure of the manifold where the target model is defined. 
In order to evaluate the universal feature quantitatively, it is necessary to 
find scaling relations of the entropy as a function of the linear system size. 
Furthermore, the scaling relations tell us how physically different systems 
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are associated with each other. When two systems have the same scaling 
relation, their eigenvalue spectra of the density matrices may be similar 
in some cases. The identification between the physically different systems 
leads to deep understanding of duality, holography, and quantum-classical 
correspondence. They are particularly important ideas in string theory and 
Quantum Monte Carlo simulation. In ergodic theory, it was a central prob- 
lem to examine what is a class of measure-preserving transformations in 
which the information entropy determines isomorphic properties. In view- 
points of symmetry and group theory, the scaling, conformal invariance, 
and hyperbolic geometry are mutually correlated. Therefore, a lot of im- 
portant concepts of physics, mathematics, and information merge together 
in the examination of the entropy. In this article, I would like to review 
some aspects related to the entropy scaling and its underlying geometrical 
structure. 



2. Entanglement Entropy and Singular Value 
Decomposition 

Let us consider a spatially <i-dimensional (dD) quantum system that is 
composed of a subsystem A and an environment B. The system is some- 
times called 'universe' or 'superblock'. The linear size of A is denoted by 
L. The entanglement entropy represents the amount of information across 
the boundary between A and B. We start with a pure state of the universe 

M=5>(*,I/)|s>li/>, (1) 

where \x) and \y) represent basis states of A and B, respectively. The density 
matrix of the subsystem A is then defined by 

PA = tr B |V) (VI , (2) 

where the symbol trs traces over degrees of freedom inside of B. Then, the 
entanglement entropy Sa is given by 

Sa = -tr A (pa In Pa)- (3) 

In order to see physical meaning of the entropy, it is better to introduce 
the singular value decomposition (SVD) of the wave function ip(x,y). The 
SVD of tf>(x, y) is defined by 



i>{x,y)=^2Ui(x)y/\lVi(y), 
i 



(4) 
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where Ui(x) and Vi(y) are the column unitary matrices, and A; is the singu- 
lar value that is positively definite. In the following, we use the normalized 
singular value pi = A;/^ ; A;. Since ^ l pi — 1, pi represents a probability 
of realization of the state labeled by I. The entanglement entropy is then 
expressed as 

S = S a = S b = -J2pi^Pi- (5) 
i 

This simple relation Sa = Sb clearly shows non-extensivity of the entan- 
glement entropy, since the volume of A is in general different from that 
of B. The object that is common in between is their boundary. Then, the 
entropy would be proportional to the boundary area. This is called area law 
scaling. In this sense, the entanglement entropy and the thermal entropy 
behave quite differently. Therefore, to confirm the area law and to find its 
violation in specific cases are two active topics. 

If we assume equal probability condition pi = 1 jm for any I , the entropy 
yields the standard Boltzmann's law 

S = - V-ln- =Inm. (6) 

' m m 
1=1 

However, the equal probability condition for the states within the subsys- 
tem A is clearly violated in low-dimensional quantum systems. It is easy 
to see this feature by, for instance, exact diagonalization calculation and 
other techniques. Thus, the difference between entanglement and thermal 
entropies is essential. Furthermore, the singular values are uniquely deter- 
mined after the decomposition, while the column unitary matrices are not. 
The universal behavior of the entropy is thus due to the presence of a set 
of the universal singular values. 

It is noted that a recent trend to examine the entanglement structure 
of quantum systems is to introduce the entanglement spectrum as well as 
the entanglement entropy. The spectrum is defined by 

Bi = -InAj. (7) 

Extensive examinations show that the spectrum is a powerful tool to char- 
acterize topological nature of the system. The nature is govened by presence 
or absense of the entanglement gap. The reason for the powerfulness comes 
from a fact that a topological phase is a phase of matter that cannot be 
described by an order parameter. 
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3. Scaling of Entanglement Entropy 

3.1. Historical Roots of Anomalous Entropy Scaling 

The thermal entropy is usually an extensive parameter. One exceptional 
case can be seen in black hole thermodynamics. 1-3 It is well known that the 
Bekenstein-Hawking entropy is proportional to the surface area A = Airr 2 
of the event horizon of a black hole: 



Here, r = 2GM/c 2 is the Schwartzshild's radius. Starting with the surface 
gravity k = GM/r 2 , the derivative form of A leads to dM — (k/8ttG)cIA. 
Due to the relativistic energy E = Mc 2 and the first law of thermodynam- 
ics dE — TdS, we obtain dS = (kc 2 /8irGT)dA. If we assume the Hawking 
temperature T = hn/2-KkBC, the Eq. (8) is derived. The information ab- 
sorbed into the black hole can not go out, if we neglect Hawking radiation. 
Then, a role of the black hole on the entropy is like enviromental degrees of 
freedom in Eq. (2), since in the definition of the entanglement entropy an 
observer in A can not access the information in B. Thus, the information 
theory based on the entanglement entropy has close connection to the black 
hole physics. 

3.2. General coodinate transformation and horizon 

Consider the metric that does not have any singularities for a static ob- 
server. However, an accerated observer may look at different spacetime 
structure. Actually, the fiat Minkowski metric can be transformed into the 
Rindler one that has event horizon. 4,5 This means that general coodinate 
transformations beyond Lorentzian ones limits spacetime region that the 
accerelated observer can access. As shown in Fig. 1, the Rindler observer is 
confined in the Rindler 'wedge' which acts as a subsystem A, while the out- 
side of the wedge represents the enviromental degrees of freedom B. Thus, 
the thermal entropy for the field propagating in this geometry behaves as 
the black hole entropy. 

Starting with the Minkowski metric (we consider ID case for simplicity) 



ds 2 = -dt 2 + dx 2 
we introduce the following transformation 



(9) 



X = ^x 2 -t 2 , 



(10) 




(11) 
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.r-»oo, x = o 

fX = const . 

1 T = const 
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X 

V Rindler wedge 

\r^-oo,x = o 



Fig. 1. Minkowski and Rindler coodinates. 

or more explicitcly 

t = XsmhT, (12) 
x = XcoshT, (13) 

where — oo < T < oo, and < X < oo. Then, the metric for the new 
coodinate is given by 

ds 2 = -X 2 dT 2 + dX 2 . (14) 

We clearly see that the accesible region is only the right Rindler wedge as 
shown in Fig. 1. 

As a tutorial example, let us consider the massless scalar field <f> propa- 
gating in the Rindler space. The action is given by 

1 = J dTdX^—g (-Ig^d^d^ . (15) 

The Lagrange equation is obtained as 

|^ = (x^ + xJL),, (16) 

and yields the following solution 

<f>(T, X) = Ae^-^ = A exp (iflXe- T ) . (17) 

By taking the inverse Fourier transform with respect to T, we see that the 
power spectrum of this field is not monochromatic 

/oo 
eff exp (inXe- T ) e lulT = (-inX) iu T(-iuj), (18) 
-OO 
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where lo 7^ 0. This leads to Planck distribution with effective inverse tem- 
perature ft — 2ir. Complete calculation of the thermal entropy, its corre- 
spondence to the Bekenstein-Hawking formula, and related aspects can be 
seen in the references. 6-10 



3.3. Area Law and its Logarithmic Violation at Criticality 

It has been extensively examined how the entanglement entropy behaves 
as functions of L and d. As already mentioned above, the most well-known 
formula is the area-law scaling 11 ' 12 

Soc(£f~\ (19) 

which tells us non-extensivity of S in contrast to the thermal entropy. This 
relation is strictly hold in gapped quantum systems. 

The violation of the area law occurs in cases of ID critical systems. 13-16 
The entropy is given by 

S=\chx-, (20) 

where c is the central charge and a is lattice cutoff. This was first derived 
from the conformal field theory (CFT), and then numerically confirmed in 
various models of statistical physics. 17-19 Away from a critical point, the 
entropy is deformed as 

S=-cA\n^, (21) 
6 a 

with correlation length £ and the number of boundary points A of A. 



3.4. Fermi Surfaces and Entropy 

The violation also occurs in models with the Fermi surfaces. 20-24 The Fermi 
surface devides the momentum space into two sectors: sets of occupied and 
unoccupied states. The electronic excitations normal to the Fermi surface 
are like those of chiral Luttinger liquid. The excitations lead to logarithmic 
correction. The entropy is then given by 

C = mo \d-i I I K • n P \ dA x dA p , (23) 
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where C is the number of excitation modes across the Fermi surface dT, 
dfl is the spatial region considered, and n x and n p are the unit normals to 
the boundaries. 



3.5. Topological Entanglement Entropy 

The physical characterization of topological orders of 2D is a hot topic in 
condensed matter physics. Usually, the topological nature can be seen as 
edge excitations and quasiparticle statistics. However, the topological order 
is manifest in the basic entanglement of the ground state wave function. 25,26 
In Z q lattice gauge theory, the entropy has a subleading term called topo- 
logical entropy in addition to the area-law scaling 

S = aL - In y/q + ■ ■ ■ . (24) 

The string-net condensed model and the Kitaev's toric code are two impor- 
tant examples of realizing the presence of the topological entropy. 

3.6. Entanglement Support of Matrix Product State 

U) 

4 "^ 

j 

Fig. 2. MPS on ID chain. Each entangled bond is characterized by a symbol fi = 
1, 2, x- The index s is local degree of freedom (s =t, X in S = 1/2 spin systems). 



In statistical physics, there is another type of entropy scaling originated 
from a recently developed variational approach to quantum systems. The 
approach optimizes the following matrix product state (MPS) ansatz 27 ' 28 

|V) = J2 tr ( A T A 2 2 ■ ■ ■ K n ) l*i*2 ■■■*„>, (25) 

where the site-dependent matrices Aj 3 have x x X dimensions, and Sj rep- 
resents local degrees of freedom (see Fig. 2). Historically, matrix product 
forms play crucial roles on exacly solvable statistical models. 

Let us consider the simplest case. We are going to expresse the singlet 
state |V>) = lt-1) ~ lit) by a particular MPS. It is clear that the direct 
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product state (local decomposition) is not exact 

10) = E tflsi)® E c 2 2 l*2> (26) 
si=T4 «2=t4 

= 44 in) + 44 in) + 44 lit) + 44 \u) ■ (27) 

Actually, we can not simultaneously take 44 = 1> c i c 2 = 0' 4 C 2 = ^, an d 
44 = — 1- However, we can introduce the following expression 

iv) - e AS1BS2 i siS2 > ( 28 ) 

S1,S2 

where the two local vectors A Sl and B S2 are taken to be 

A^=(xy),A^=(zw), (29) 

B^^(y ),B± = ^( W ), (30) 

xw — yz \ —x J xw — yz \ —z J 

for scalar variables x, y, z, and w. This expression means that the matrix 
dimension \ provides quantum correlation between neighboring spins. 

The MPS state is known to be exact for gapped ID systems such as 
the valence bond solid (VBS) state. Furthermore, the MPS is foundation 
of the density matrix renormalization group (DMRG) method that is a 
powerful numerical technique for ID quantum systems. 29,30 DMRG devides 
the system into two parts, and then truncates the density matrix eigenstates 
with the small eigenvalues. Therefore, DMRG is based on entanglement 
control. The restriction of the matrix dimension upto \ means that the 
amount of information described by this MPS is given by 

x 

^ x = -E A ' lnA '' ( 31 ) 
i=i 

where Ai > A2 > • • • > A x . If we assue that A* J is a scalar variable (x = 1), 
Eq (25) gives a local approximation as I have already mentioned. Taking 
a sufficiently large x value gives us an asymptotically exact wave function. 
Thus, the matrix dimension \ represents how precisely we can take long- 
range quantum correlation. When we apply this variational method to ID 
critical systems, the half-chain entropy (.4 = 1) behaves as 31-33 

ck . 
~6~ 

where the exponent k is defined by 

6 



5 x = -ln X , (32) 



12/c+ 1 



(33) 
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Comparing Eq. (32) with Eq. (21), we know that 

£ = X K , (34) 

which means that \ is actually related to the length scale £. For instance, 
in the Ising universality class c = 1/2 

k ~ 2, (35) 

S x ~\\^x- (36) 

I have found an old DMRG paper that addresses the correlation length £ 
as a function of the truncation number. 34 The result for ID free fermion 
with the central charge c = 1 is given by 

1 1 



s lnll-fcx^l k X ' (37) 

with /? <~ 1.3 and k <~ 0.45. This is really consistent with the above argument 
where k <~ 1.33 for c = 1. 



4. Holographic Entanglement Entropy: Connection 
between Scaling and Hyperbolic Geometry 

4.1. Anti-de Sitter space 

In superstring theory, it has been an interesting topic to find the fundamen- 
tal mechanism of the anti-dc Sitter space (AdS) / CFT correspondence. 
This is holographic correspondence between a quantum theory, CFTd+i, 
and general relativity on the AdS space, AdSd+2 (Here, CFT^+i is ab- 
breviation of spatially c?-dimensional conformal field theory with one time 
axis). The AdS space is (d + 2)-dimensional hyperbolic surface embedded 
in one higher dimensional flat space E d+1 ' 2 

-X 2 _ 2 -Xl 1+ J2x 2 k = -l 2 , (38) 

k=0 

(see Fig. 3). This space has two time- like directions A_ 2 and A_i. The 
metric of the AdS space can be represented as 

ds 2 = g^dx^dx" (39) 

= —(dz 2 + r ]ij dx l dx'), (40) 

where I is the curvature of this space. The index z is called radial axis, and 
i runs over 0, 1, 2, d. The AdS space is a vacuum solution of the Einstein 
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equation 



- \g^R + .9^A = kT"" = 0, (41) 



with negative cosmological constant A = —d(d + l)/2l 2 . 

4.2. Killing equation and conformal invariance 

The Killing equation on the bulk AdS space approaches the conformal one 
at the boundary z — > 0. Let us start from infinitesimal transformation 

^ =x' t +^(x), (42) 

and then the new line element is given by 

ds' 2 = ds 2 + (V„&, + V„£„) da^dx", (43) 

where the covariant derivative is defined by 

V„& - 0„£„ - r*„&, (44) 

with the Christoffcl symbol 

= \<T (d»g„ r + d v9llT ~ d T9ltv ) . (45) 

Here, the isometory ds' 2 = ds 2 yields the Killing equation V M £„+V„£ M = 0. 
We consider the following transformation 

C(z,x)^e i (x), (46) 
C(z,x)^zC(x), (47) 

where the boundary does not move after the transformation. Substituting 
these equations and the AdS metric into the Killing equation, we actually 
obtain the conformal Killing equation at the boundary 

di€j + djei = CVij- (48) 

Thus, the CFT lives on the boundary of the AdS space. Actually, it is 
possible for various models to prove 



S 5 (I 

ex P 



6<f)(x) 6<j>(x') ^\ 2k 2 



<* TZ ZTTa"' (49) 



b=0 



\X — X 



with the scaling dimension A of the CFT and / denoting the classi- 
cal Einstein-Hilbcrt action. This is so called Gubser-Krevanov-Polyakov- 
Wittcn relation. 35-38 
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4.3. Geometric interpretation of entanglement entropy 

An important guiding result to understand the physical background of the 
AdS/CFT correspondence is the so called Ryu-Takayanagi's formula for the 
entanglement entropy 39-41 



Areafa) 
4G : 



(50) 



where ja is the minimal surface whose boundary is given by the manifold 
d-fA — dA, and G is the Newton constant of gravity in the AdS space. This 
is general extension of the Bekenstein-Hawking entropy in Eq. (8). 




(a) 




Fig. 3. Two representations of 2D hyperbolic space and geodesic line: (a) Poincarc disk 
model, (b) discrete model. In (b), Dxy may not look like a geodesic line, but this is 
due to ambiguity of lattice discretization. The original AdS metric in Eq. (40) can be 
transformed into ds 2 = (dr log2) 2 + (2~ T / l dx) 2 . This form is comparable to Fig. (b). 



For d = 1, Area(7yi) is the geodesic distance Dxy between X = 
(— L/2,a) and Y = (L/2,a) for a — > as shown in Fig. 3. These points 
are located on the boundary of the AdS space. Let us derive an explicit for- 
mula of Dxy, since this quantity plays a crucial role on the later discussion. 
The geodesic equation is given by 

' +r%^-^ = o, (5i) 



dr 2 ^ dr dr 
and the solution is the half-cycle 



(xV 2 ) = (x,z) = | (cos 0, sin 0), (52) 
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where d6/dt = sin 6*, e < 9 < ir — e, and (L/2)sine = a. By using the 
equation for the geodesic line, we calculate Dxy as 

Dxy = 2 C' 2 -d6sj \d e zf + (d s x) 2 



/In 



1 + cos e 
1 — cos e 



By combining this with the Brown-Henneaux central charge 42 

°"55' <54 » 

we can derive the CFT result S — (c/3) h\(L/a) for L ^> a. This calculation 
can be extended to higher dimensional cases. 

It is noted that physical meanings of curved spaces are also explored 
in statistical mechanics. Hyperbolic deformation of interaction and sine- 
squared deformation are two major approaches. 43-46 

5. Tensor Networks and Extra Dimension 
5.1. From Matrix Product to Tensor Product 




Fig. 4. Tensor network structure on 2D square lattice. The tensor T^ vaT with linear 
dimension \ is defined on each lattice site. The quantum entanglement between neigh- 
boring sites is represented by the indices fi, v, a, and r. 
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As we have already mentioned, the matrix dimension \ of the MPS 
state represents how strong quantum entanglement appears. However, the 
MPS state is appropriate only in ID gapped cases. Direct extention of this 
approach to higher dimensions is to construct tensor networks (or some- 
times called tensor product state, TPS, and projected entangled-pair state, 
PEPS). 47-52 We show a schematic viewgraph of TPS in Fig. 4. We define a 
tensor, T 1 ^..., on eacn s ^ e 3- The num ber of the indices, /j,, u, A, • • • , of each 
tensor corresponds to that of connected bonds. If each bond is maximally 
entangled, the state of the bond can be represented by a and its ancila a 



Then, the entropy for each bond is given by Sbond = —tr (phip) = logx, 
where the trace is taken over the ancila degree of freedom. The entanglement 
entropy of the subsystem marked by a dotted line in Fig. 4 is given by 



where N is the number of the entangled bonds, L is linear size, and a is 
lattice constant. This relation is nothing but the area law scaling. Thus, 
we know that the TPS state is suitable for (i-dimensional gapped cases. On 
ther other hand, x ~ L in critical cases, and further improvement of the 
network structure may be required. This is the purpose of the next section 
where tree tensor and MERA networks are introduced. 

5.2. Tree Tensor Networks and Multiscale Entanglement 

Renormalization Ansatz: Hierarchical Tensor Network 

Constructing efficient variational optimization methods of MPS and TPS 
is of practical importance in condenced matter physics. An advantage of 
these methods is to keep the area-law scaling of the entropy, and thus is 
suitable to gapped (i-dimensional systems. The application of tensor net- 
works to critical systems leads to the tree tensor network (TTN) and the 
multiscale entanglement renormalization ansatz (MERA). 53-60 They are 
spatially (d+ l)-dimensional hierarchical networks that are also compatible 
with real-space renormalization group. In particular, the MERA network 
forms discrete AdS space as shown in Figs. 5 and 6. 61 ' 62 Although strict 
correspondence between AdS/CFT duality and MERA is not constructed 
in the present stage, they are very similar. When we express a critical sys- 
tem by using TPS, we need large tensor dimension \ ~ L. On the othere 




(55) 




(56) 
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X Y 

Fig. 5. Schematic MERA network. Triangles and diamonds arc isometory tensors and 
discntanglcr tensors, respectively. Each bond connecting a triangle with a diamond rep- 
resents tensor product. For comparison, Fig. 3 is combined with the network. 




Fig. 6. Deformation of MERA network in Fig. 5. 



hand, each tensor in the hierarchical network of MERA has relatively small 
dimension. Thus, we may say that MERA constructs more classical-like 
states in the holographic space. 
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6. Compactification and Entropy 

6.1. Bulk/edge Correspondence and Compactification 




Fig. 7. Compactification of open 2D sheet into ID system 

AdS/CFT is a kind of bulk/edge correspondence on non-compact man- 
ifold. Compactification of non-compact space is an alternative way to find 
correspondence between physically different systems. The comparison be- 
tween these ways may be useful for further examination of the entangle- 
ment entropy. The concept of the compactification is schematically shown 
in Fig. 7 where an open 2D sheet is rolled up. When our spatial resolu- 
tion is worse than the compactification radius, this is effectively a ID sys- 
tem. Unfortunately, we can not directly define the entanglement entropy 
in the classical space, since entanglement represents quantum correlation. 
However, according to the Ryu-Takayanagi's formula, a geometric quantity 
would correspond to the entropy. Then, the problem is to find geometrical 
meaning of the entanglement entropy in the non-compact manifold before 
compactification. In the following subsections, some historical review and 
my recent work closely related to this aspect are presented. 

6.2. Field theory with Extra Dimension 

Let us look at compactification of a massless scalar field model. The met- 
ric is given by the flat 4D Minkowski one (i = 0,1,2,3) plus additional 
component dx 4 = ad9 

*— (o"0- (57) 

where the radial direction is represented by 9, and a is its radius. We assume 
a cylindrical boundary condition for the 8 direction, and the others are not. 
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Then, we have < 9 < 2ir. The scalar field (j>(x, 9) is discretized as 

<t>(x,0) = -±=Y,^{x)e ml \ 
V 27r 



(58) 



where </>_„ — 4>n- Substituting it to the action, the action after integrating 
over 9 degree of freedom is given by 

I= — I d d x I d8^—g { -^d a (t)d u 4> 1 (59) 



a J d d a 



n>l 



■ (60) 



Therefore, in addition to the zero mass mode cfio, we obtain the massive 
modes <j) n with mass M n = n/a. They are called Kaluza-Klein modes. In 
the classical field theory which aims dimensional reduction, a should be 
taken to be small enough. However, in the following, we are interested in 
cases of various a values. For large a values, these modes tend to be gaplcss. 
Thus, this indicates that the enough internal degree of freedom is necessary 
to describe low-energy excitations. This feature can also be seen in the 
tensor network formulation of quantum states. 



6.3. VBS/CFT correspondence 

Recently, it has been shown that the entanglement spectra of the VBS state 
on 2D lattices are closely related to a thermal density matrix of a holo- 
graphic spin chain whose spectrum is reminiscent of that of the spin- 1/2 
Heisenberg chain. 63,64 Here, we briefly touch on an idea of this correspon- 
dence based on the tensor product. It might be possible to have a different 
view based on open/close string duality. We start with the 2D VBS state 
on M x N lattice (vertical x horizontal) expressed by 

where /„ ~ («i,n, *2,n, *M.n)- Let us consider a cylindrical boundary con- 
dition that the vertical direction is rolled up, while the horizontal axis 
remains open. Then, we can introduce the coefficient a defined by 

°I = Yl ^Ai^AiAs ' ' ' ^A^UiV-i^Aw-i' ( 62 ) 
A 

where A„ = {ai tTl , «2, n , aiM.n) and aj^ n = 1, 2, \ with bond dimension 
X of the original VBS. The periodic boundary condition along the vertical 
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direction due to the cylinder form is thus expressed by 

M 

B t-^n= tr Tl^-^ (63) 
j = l 

and its boundary terms L 1 ^ and J2^_ . Clearly, Eq. (62) is a MPS form. 
When we devide this system into two parts, the state can be expressexd 
by the Schmidt decomposition 

where L Ia = L h B h ■ ■ ■ B h and R Ib = B Il+1 ■ ■ ■ B In ~ 1 R In . The dimension 
of A is x M ■ Originally, the VBS state has relatively small dimension \ 
due to the presence of the Haldane gap. However, x M becomes very large, 
leading to critical behavior. I think it's interesting to compare this idea 
with the compactification discussed in the previous subsection. 



6.4. Suzuki- Trotter Decomposition 

The readers expert for quantum Monte Carlo (QMC) simulation may be 
aware of the similarity between Suzuki- Trotter decomposition (STD) and 
compactification. 65 The correspondence between classical and quantum sys- 
tems are induced by the following STD 

e x+Y = lim (eie» 

M— >oc V 

for non-commutative operators X and Y. In some cases, the right hand side 
can be treated exactly. 

It is well-known that the partition function of the transverse-field Ising 
chain can be mapped onto that of the 2D classical anisotropic Ising model. 
The transverse-field Ising Hamiltonian is given by 

L L 

i=l i=l 



(65) 



Let us derive the classical model with use of the STD. We introduce the 
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partition function Z — Tre @ H , and decompose it into M slices as 



M 



exp ( ——H I cxp ( ——H' 



(3 



M 



M 



IW> 



J2 ft (Ml <*P (~Ji H o) ^P (-^^') IW+ 

{<Tl} {<TM}fe=l V / \ / 



(67) 

i}> (68) 



M 



- e n-P ^e^i 



{cti},...,{ctm } fe=l 



1=1 



/3 



X (K}|CX P ( -^H'J \{cTk+l}) 



(69) 



Here, we have imposed the cylindrical boundary condition, {<Ji} = {<xm+i}- 
Thus, we may say that the mathematical processes of the STD zoom up 
internal degrees of freedom on each site, which is analogous to the previous 
sections. For a, a' = ±1, we have a relation 



{<j\e^ Xa " \a') = A exp {-^' In (^tanh-^A^ J 
with A = v /(l/2)sinh(2/3/M)A. Thus, Z is expressed as 

fc=l I i=l J 

xexpj-iXJa^ln^tanhAA^J 
and finally the effective Hamiltonian is given by 

L M 

#eff = e E ( j i^^+i + , 

j=l k=l 

where Ji = -J/M and J 2 = -(1/2/3) In (tanh(£A/M)). 



(70) 



(71) 



(72) 



6.5. Entropy scaling, quantum- classical correspondence, 
and hyperbolic geometry hidden in image processing 
based on SVD 

According to the previous discussion, let us finally examine what is a geo- 
metrical object corresponding to the entropy on the classical spin systems 
before compactification. That is a snapshot of a particular spin configura- 
tion at criticality. 
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We regard the orignal image data of the snapshot with M x N pixels 
as a matrix ip(x,y) (1 < x < M, 1 < y < N). We assume that i/j is real. 
Then, it is possible to introduce the density matrices 

px(x,x') = ^2ip(x,y)il>(x',y), (73) 
v 

py{y, y') = J2 f M 1 ' f ')■ ( 74 ) 

X 

Their L non-zero eigenvalues are the same (L = min(M, N)). According to 
Eqs. (4) and (5), we can define the entropy. In order to intuitively under- 
stand scaling relations which this snapshot entropy satisfies, let us imagine 
a snapshot of the classical 2D isotropic Ising spin system at criticality. 
There, the spin configuration is fractal, since various sizes of the ordered 
clusters coexist due to the critical fructuation. This also means that vari- 
ous length scales are mixed. We can pick up many fractal subsystems from 
the original snapshot, but their patterns themselves arc different with each 
other. Therefore, after correcting all of the possible patterns, they would 
cover the total information of the partition function. The information of 
thermal fructuation contained in the partition function is represented by 
the differen patterns. In that sense, only one snapshot is necessary at crit- 
icality. According to the STD, we expect that there exists a quantum ID 
system that is transformed into 2D isotropic classical Ising model. There- 
fore, the snapshot entropy should obey the scaling relation equal to that 
of the critical ID quantum systems. Actually, I have recently shown that 
this conjecture is at least correct in the 2D isotropic classical Ising model 
on L x L lattice. 66 I have obtained the following results at T c 

lnL-2, (75) 

+ (76) 

for sufficiently large L and a positive 7 value. They are comparable to CFT 
and MPS results. Furthermore, when we visualize each layer of the SVD 

^ l \x,y) = U l {x)^K l V l {y), (77) 

we know that ip( l \x,y) has its own length scale. The images of tp^(x,y) 
with larger singular values contain more global spin structures. Extensive 
examinations suggest that the index I corresponds to the radial direction 
of the discrete AdS space. The result means that the coarse graining of 



1=1 

S x = -^A ; lnA* 



1=1 
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the image by tuncating small singular values correspondes to a flow from 
the boundary to the bulk on the AdS space. This would be related to the 
holographic renormalization group. 

The entropy of more realistic images also obey a clear scaling relation 
analogous to the entanglement support of MPS. 67 This is surprizing, since 
we think it's impossible to determine 'the central charge' of the realistic 
images. We need more detailed analysis of this type of quantum-classical 
correspondence. 

7. Summary 

In this article, I have reviewed various scaling relations of the entanglement 
entropy and related topics. Since the entropy represents universality of the 
model considered, we can obtain more global viewpoints beyond the model 
itself. In particular, the scaling is a very powerful tool to look at quantum- 
classical correspondence. Finally, we should be carefull to a fact that the 
quality of the information is given by the wave function, not the entropy. 
In this respect, we need comprehensive study of the scaling analysis and 
numerical techniques of optimizing tensor networks. 
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